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Abstract. We present a new approach to the minimum-cost integral 
flow problem for small values of the flow. It reduces the problem to the 
tests of simple multi- variable polynomials over a finite field of characteris- 
tic two for non-identity with zero. In effect, we show that a minimum-cost 
flow of value k in a network with n vertices, a sink and a source, integral 
edge capacities and positive integral edge costs polynomially bounded 
in n can be found by a randomized PRAM, with errors of exponentially 
small probability in n, running in 0(fclog(fcn) + log 2 (fcn)) time and using 
2 fe (fcn)°( 1 ) processors. Thus, in particular, for the minimum-cost flow of 
value O(logn), we obtain an RNC 2 algorithm. 

1 Introduction 

The maximum network flow problem is a well known fundamental problem in 
algorithms and optimization with plenty of important applications [1,7,8,16]. 
It is known to be P-complete even in its integral version provided that the edge 
capacities are exponentially large in the size of the network [13]. The minimum- 
cost flow problem is a well known important generalization of the maximum flow 
problem [1, 8, 10, 16]. The objective is to compute a maximum flow of minimum 
cost in a directed graph where each edge is assigned a cost. For a flow / in a 
directed graph (V, E), the cost of / is simply J2 e eE f(e)cost(e). 

The prospects for designing a fast and processor efficient parallel algorithm, 
in particular an NC algorithm [19], for maximum integral flow or minimum-cost 
integral flow are small. The fastest known parallel implementations of general 
maximum flow and/or minimum-cost flow algorithms achieve solely a moderate 
speed up and still run in f2(n a ) time, where a is a positive constant, see [2, 12]. 

The situation changes when the edge capacities or the supply of flow as well 
as edge costs are substantially bounded. For example, if the edge capacities 
and edge costs are bounded by a polynomials in n, both problems admit RNC 
algorithms. Then, the maximum integer flow problem admits even an RNC 2 
algorithm [17, 18, 20] while the minimum-cost integer flow problem admits an 
RNC 3 algorithm [18]. At the heart of the aforementioned RNC solutions is 
the randomized method of detecting a perfect matching by randomly testing 
Edmonds' multi- variable polynomials for non-identity with zero [6, 14, 17, 20]. 



When the flow supply is relatively small, e.g., logarithmic in the size of the 
network or a poly-logarithmic one, then just an NC implementation of the basic 
phase in the standard Ford-Fulkerson method [1, 7-9, 16] yields an NC algorithm 
(NC 3 when the supply is logarithmic) for maximum integer flow that can be 
extended to an NC algorithm for minimum-cost integer flow (when edge costs 
are polynomially bounded). The number of processors used corresponds to that 
required by a shortest path computation. 

In this paper, we present a new approach to the minimum-cost integral flow 
problem for a small value k of the flow. We directly associate a simple polynomial 
over a finite field with the corresponding problem of the existence of k mutually 
vertex disjoint paths of bounded total length, connecting two sets of k terminals 
in a directed graph. By using the idea of monomial cancellation, the latter prob- 
lem reduces to testing the polynomial over a finite field of characteristic two for 
non-identity with zero. We combine the DeMillo-Lipton-Schwartz-Zippel lemma 
[5, 21] on probabilistic verification of polynomial identities with parallel dynamic 
programming to perform the test efficiently in parallel. Additionally, we use the 
isolation lemma to construct the minimum-cost flow [17,20]. 

In effect, we infer that a minimum-cost flow of value A: in a network with n 
vertices, a sink and a source, integral edge capacities and positive integral edge 
costs polynomially bounded in n can be found by a randomized PRAM, with 
errors of exponentially small probability in n, running in 0(k log(kn) +log 2 (fcn)) 
time and using 2 k (kn)°^ processors. Thus, in particular, for the minimum-cost 
flow of value O(logn), we obtain an RNC 2 algorithm. 

Related work. For the RNC algorithms for the related problem of minimum- 
cost perfect matching see [14, 17, 20]. For the comparison of time and substantial 
processor complexities of prior RNC algorithms for the minimum-cost flow see 
page 7 in [18]. The fastest of the reported algorithms is not an RNC 2 one even 
when the flow supply and thus the edge capacities are logarithmic in the size of 
the network. The idea of associating a polynomial over a finite field to the sought 
structure has been already used by Edmonds to detect matching [6] and then in 
several papers presenting RNC algorithms for perfect matching construction [14, 
17, 20]. It appears in several recent papers that also exploit the idea of monomial 
cancellation [3,4,15,22]. 

Organization. In the next section, we comment briefly on the basic notation 
and the model of parallel computation used in the paper. In Section 3, we derive 
our fast randomized parallel method for detecting the existence of k mutually 
vertex disjoint paths of bounded total length connecting two sets of k terminals 
in a directed graph. In Section 4, we generalize the method to include edge 
costs which enables us to replace the total length bound with the total cost one. 
In section 5, we show a straightforward reduction of the minimum-cost integer 
flow problem parametrized by the flow value to the corresponding disjoint paths 
problem which enables us to derive our main result on detecting minimum-cost 
small flows in parallel. 



2 Terminology 



For a natural number n, we let [n] denote the set of natural numbers in the 
interval [1, n]. The cardinality of a set A will be denoted by \A\. 

We assume the standard definitions of flow and flow value in a network 
(directed graph) with integral edge capacities, a distinguished source vertex s 
and a distinguished sink vertex t (e.g., see [7]) . 

For the definitions of parallel random access machines (PRAM), the classes 
NC and RNC and the corresponding notions of NC and RNC algorithms, the 
reader is referred to [19]. 

The characteristic of a ring or a field is the minimum number of 1 in a sum 
that yields 0. A finite field with q elements is often denoted by F q . 

3 Connecting vertex-disjoint paths 

It is well known that the maximum integral network flow problem with bounded 
edge capacities corresponds to a disjoint path problem (cf. [7]). In Section 5, we 
provide an efficient parallel reduction of the minimum-cost integral flow problem 
parametrized by the flow value to a parametrized disjoint path problem. This 
section is devoted to a derivation of a fast randomized parallel method for the 
decision version of the parametrized path problem. 

Let L = (V, E) be a network in a form of a directed graph with n vertices, 
among them a distinguished set X = {x\, Xk} of k source vertices and a 
disjoint distinguished set Y = {j/i, j/fe} of k sink vertices. 

A walk in L is a sequence of vertices V\, V2, Vi of L such that for j = 
1, I — 1, (vj,Vj + \) G E, v\ is in X, v 2 , arc in V"\ (XUY), vi is in Y. The 
length of the walk is I — 1. In other words, a walk is just a (not necessarily simple) 
path starting from a vertex in X, having intermediate vertices in V \ (X U Y), 
and ending at a vertex in Y. 

A proper set S of walks in L is a set W of k walks of total length < k(n — 1), 
each with a distinct start vertex in X and a distinct end vertex in Y. 

A signature of a proper set S of walks is the pair that is smallest in 

lexicographic order such that the two walks that start at Xi and Xj respectively 
intersect, and the first intersection vertex of these two walks is the first intersec- 
tion vertex of the walk starting from a;, with any walk in S. 

Note that walks in S are pairwise vertex disjoint iff the signature of S is not 
defined. 

We define the transformation <p on S as follows. If S has the signature 
then (f> switches the suffix of the walk starting at xi with that of the walk starting 
at Xj at the first intersection vertex of these two walks. See Fig. 1. Otherwise, 
if the signature of S is not defined then <j> is an identity on S. 

Observe that if the signature of S is defined then 4>(S) has the same signature 
as S and (j>(S) ^ S. The first observation is immediate. To show the second one 
it is sufficient to note that <fi(S) = S holds iff <f> transforms the two walks which 
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Fig. 1. An example of a proper set S of walks and the companion proper set 0(5*) of 
walks. 

yield the signature of S onto themselves. The latter is however impossible since 
they have different start vertices and different end vertices. Note also that the 
walks in <j>(S) have the same total length as those in S. 

It follows that <f) is an involution on sets of proper walks of total length I, 
i.e., 4>{4>{S)) = S holds for any proper set S of walks of total length I. 

For the network L and I £ [k(n— 1)], let Fl^i be the family of all proper sets of 
k walks of total length < / in L. Assign a distinct variable x e to each edge e in L. 
For a walk W € F^j, let be the monomial, where x e has multiplicity equal 
to the number of occurrences of e in W. Next, let Ql.i denote the polynomial 

Lemma 1. For the network L and I G [k(n — 1)], there is a proper set of k 
mutually vertex- disjoint walks of total length < I in L iff Ql.i is not identical to 
zero over a field of characteristic two. 

Proof. Fl,i can be partitioned into the family F\ ; of sets S of walks such that 
(j>(S) = S and the family F 2 L l of sets 5 of walks such that 0(5) ^ S. The 
polynomial X^seF 2 Ilwes is identical to zero over a field of characteris- 
tic two since for each S £ F\ the monomials Jlvves an d Iliye</>(S) 
contain equal multiplicities of the same variables and (j)(<f)(S)) — S so S and 
<fi(S) can be paired. On the other hand, since each set S of walks in F^ ( consist 
of mutually vertex-disjoint walks, the monomials Jlwes m the polynomial 
XseFi , IlwGS M\y are in onc-to one correspondence with S and thus arc unique 
provided that the walks in S are simple paths. Now, it is sufficient to observe that 
mutually vertex-disjoint walks can be always trivially pruned to corresponding 
mutually vertex-disjoint simple paths. 

□ 



To warm up, we prove the following lemma on sequential evaluation of Ql- 



Lemma 2. Ql : i can be evaluated for a given assignment of values over a field 
F 2 o(io g n) of characteristic two in 0(k 3 n A + 2 k k 4 n 3 ) time. 

Proof. For B C Y, I £ [(n - 1)\B\], we consider the family W t (B) of all sets S 
consisting of \B\ walks connecting \B\ distinct sources in {a?i, aJisi } with the 
\B\ distinct sinks in B so that the total length of the walks is exactly I. Next, 
we define the polynomial Q P (B) as J2sew,(B) Ylwes M\y- Note that Ql,i = 

H P =k Qp(Y) an d l < nk. 

On the other hand, for p £ [k(n — 1)], x £ X and z £ V \ X, we consider the 
set W p (x, z) of walks of length p in L that start at x and end at z. Let Q p (x, z) 
be the polynomial J2wew p (x,z) M w- 

We have the following recurrence for a nonempty subset B of Y and p £ 
[\B\(n-l)}: 

Qp(B) = Y, E Q P - q (B\{y})Q g (x lBl ,y). 

VGB q e[\B\(n-l)-\B\ + l] 

Next, we have also the following recurrence for x £ X, z £ V \ X, and 
q £ [k(n - 1)] : 

Q q (x,z)= Qq-l{x,U)X {UjZ) . 
ueV\(XUY)&z(u,z)eE 

We have also Qi(x, z) = £( x ,z) if ( x i z ) € E, and otherwise Q\(x, z) — 0. Con- 
sequently, we can evaluate all the polynomials Q q (x, z) by the second recurrence 
in 0(k 2 n 3 ) time. 

Now, by using the first recurrence and setting Q p (fb) to 1 in the field, we can 
evaluate all the polynomials Q P (B) in the increasing order of the cardinalities 
of B in 0(2 k k 3 n 2 ) time. 

By Ql.i — J2 P =k Qp(Y) an d I < nk, we conclude that Ql,i can be evaluated 
in 0{k 3 n A + 2 k k i n 3 ) time. 

□ 

We can partially parallelize the sequential evaluation of Ql,i in order to 
obtain the following lemma. 

Lemma 3. Ql,i can be evaluated for a given assignment of values over a field 
F 2 o(ia Sn ) of characteristic two in 0(klogn + log 2 n) time by a CREW PRAM 
using 0(k 2 n 5 + 2 k k 4 n 3 ) processors. 

Proof. We generalize the definition of the set W q (x,z) and the corresponding 
polynomial Q q (x,z) to include arbitrary start vertex x £ V \ Y, requiring z £ 
V \ X as previously. Then, we can evaluate Q q {x,z) — ^2\yeWi(x z) Mw for 
x £ V \Y and z £ V\X, for q £ [k(n — 1)] by the following standard doubling 
recurrence for q > 2 : 

Q q {x,z)= Y Q[q/2]{x,y)Qi q /2}(y,z) 

y ev\{xuY) 



At the bottom of the recursion, we have Q\(x,z) — X( x>z ) if (x, z) e E, 
otherwise Qi(x, z) = 0. It follows that all Q q (x,z) for x € V\Y, z e V\X, and 
g G [fc(n — 1)] can be evaluated in a bottom-up manner in 0(log 2 n) time by a 
CREW PRAM using 0{kn 4 ) processors. 

Recall the first recurrence from the proof of Lemma 2. When the polynomials 
Q g (xi,yj) for Xi <G X and yi € Y are evaluated, we can evaluate in turn the 
polynomials Q P (B), where B C Y, p € [|i?|(n — 1)] in A; phases in the increasing 
order of the cardinalities of B by this recurrence. It can be done in 0(fc(log k + 
logn)) = O(fclogn) time by a CREW PRAM using 2 k k 3 n 2 processors. 

By Ql,i = ^2 p —kQp0^)i ' — an d k < n, we conclude that Ql,i can be 
evaluated in 0(fc logn + log 2 n) time by a CREW PRAM using 0(k 2 n 5 + 2 k k 4 n 3 ) 
processors. 

□ 

The following lemma on polynomial identities verification has been shown 
independently by DcMillo and Lipton, Schwartz, and Zippcl. 

Lemma 4. [5, 21] Let Q(xi,x 2 , x m ) be a nonzero polynomial of degree d over 
a field of size r. Then, for /i, fi, ■■■,f m chosen independently and uniformly at 
random from the field, the probability that Q(fi, /2, f m ) is n °t equal to zero 
is at least 1 — ^ . 

Note that the polynomial Ql,i is of degree I not larger than k(n— 1) < n 2 . We 
can use Lemma 4 with a field F 2 a io g n of characteristic two to obtain a randomized 
test of the polynomial Ql,i for not being identical to zero with one side errors. 
For sufficiently large constant c, the one side errors are of probability not larger 
than a constant smaller than 1. By performing 0(n) such independent tests, the 
probability of one side errors can be decreased to exponentially small in n one. 

By Lemma 3, the series of the tests can be performed in 0(fclogn + log 2 n) 
time by a PRAM using 0(k 2 n 6 + 2 k k 4 n 4 ) processors. By Lemma 1, these tests 
verify if there is a proper set of mutually vertex-disjoint walks of total length 
< I in the network L. The latter in turn is equivalent to the existence of k 
mutually vertex-disjoint paths of total length < I connecting X with Y in L by 
the definition of a proper set of walks in L. Hence, observing that each walk can 
be trivially pruned to a simple directed path with the same endpoints, we obtain 
our main result. 

Theorem 1. The problem of whether or not there is a set of k mutually vertex- 
disjoint simple directed paths of total length < I connecting X with Y in the 
network L can be decided by a randomized CREW PRAM, with one-sided errors 
of exponentially small probability in n, running in 0(k\ogn + log 2 n) time and 
using 0(k 2 n e + 2 k k 4 n 4 ) processors. 

4 Vertex-disjoint connecting paths of bounded cost 

In this section, we shall consider a more general situation where there are a 
positive integer C and a cost function c assigning to each of the m edges e in 



the network L a cost c(e) G [C] . The cost of a walk or a path is simply the sum 
of the costs of the edges forming it (the cost of an edge is counted the number 
of times it appears on the walk or path). We would like to detect a proper set 
of k walks in L that achieves the minimum cost. 

For this reason, we consider the following generalization of the polynomial 
Ql,i- For U e [mC], let Hl,u be the set of all proper sets of walks in the 
edge-costed network L that have total cost not greater than U. Next, for a walk 
W in L, as previously, let M w be the monomial which is the product of x e 
over the occurrences of edges e on W. The polynomial CQl,u is defined by 

J2seH LiU Uwes M w- 

By using the proof method of Lemma 1, we obtain the following counterpart 
of this lemma for CQl,u- 

Lemma 5. For the edge-costed network L, there is a proper set of k mutually 
vertex- disjoint walks of total cost < U in L iff CQl,u is n °t identical to zero 
over a field of characteristic two. 

Next, we obtain the following counterpart of Lemma 3 for CQl,u- 

Lemma 6. CQl.u can be evaluated for a given assignment f of values over a 
field F 2 o(io gre ) of characteristic two in 0(k log(Cn) + log 2 (Cn)) time by a PRAM 
using 0(k 2 C 5 n 10 + 2 k k 4 C 3 n e ) processors. 

Proof. The proof reduces to that of Lemma 3. We replace each directed edge e 
of cost c(e) € [C] in the network L by a directed path of length c(e) introducing 
c(l) — 1 additional vertices. With each edge on such a path, we associate a 
variable. We assign f{x e ) to the variable associated with the first edge on the 
path replacing e, and just 1 of the field to the variables associated with the 
remaining edges on the path. 

The resulting network L' is of size 0(Cn 2 ). Let Hl',u be the family of all 
proper sets of k walks of total cost < U in the network V . We can evaluate 
the polynomial Ql',u — ^2seH L , v Y\wes m parallel analogously as Ql,i in 
the proof of Lemma 3. It remains to observe that the value of CQl.u under the 
assignment / is equal to that of Ql',u under the aforementioned assignment. □ 

Now, we are ready to derive our main result in this section. 

Theorem 2. The minimum cost of a set of k mutually vertex- disjoint simple 
directed paths connecting X with Y in the network L with edge costs in [C] can 
be computed by a randomized CREW PRAM, with errors of exponentially small 
probability in n, running in 0(k log(Cn)+log 2 (Cn)) time and using 0(k 2 C 6 n 13 + 
2 k k 4 C 4 n 9 ) processors. 

Proof. The minimum cost of the sought set of vertex-disjoint paths is in [Cn 2 ]. 
Hence, by Lemma 5, it is sufficient to test the polynomials CQl.u for non- 
identity with zero for all U € [Cn 2 ] in parallel. By applying Lemmata 4 and 6 in 
a manner analogous to the proof of Theorem 1, we conclude that it can be done 
by a randomized CREW PRAM, with one-sided errors of exponentially small in 
n probability, running in 0(k\og(Cn) + log 2 (Cn)) time and using 
0(Cn 2 xnx (k 2 C 5 n w + 2 fe fc 4 C 3 n 6 )) processors. □ 



5 Finding vertex-disjoint connecting paths 



A straightforward approach of extending our randomized parallel method for 
deciding if there is a proper set of k mutually vertex-disjoint walks (of a bounded 
total cost) between two sets of vertices of cardinality fc to include the finding 
variant could be roughly as follows. In parallel, for each fc-tuple of respective 
neighbors of the k start vertices in X, replace the set of start vertices by the 
fc-tuple and apply our method recursively to the resulting network. If the test is 
positive, the first edges on the walks are known, and we can iterate the method. 
The problem with this approach is that its recursive depth is proportional to the 
maximum length of a walk in the resulting set of mutually vertex-disjoint walks 
between X and Y. 

Also, it is not clear how one could implement a straightforward divide-and- 
conquer approach of guessing intermediate vertices in order to find a set of fc 
mutually-vertex disjoint walks of a given cost efficiently in parallel. 

We need more advanced methods to obtain a very fast parallelization of the 
finding variant. We shall modify the edge cost in the network L in order to use 
the so called isolation lemma in a manner analogous to the RNC method of 
finding a perfect matching given in [20] . 

Lemma 7. (The isolation lemma [20]). Let F be a family of subsets of a set 
with q elements and let r be a non-negative integer. Suppose that each element s 
of the set is independently assigned a weight w(s) uniformly at random from [r], 
and the weight of a subset S in F is defined as w(S) = J2xgs w ( x )- Then, the 
probability that there is a unique set in F of minimum weight is at least 1 — * . 

Corollary 1. For each of the m edges e in the network L, modify its cost c(e) to 
c'(e) = c{e)rm + w{e), where the weight w(e) is drawn uniformly at random from 
[r]. Then, the probability that there is a unique minimum-cost set of mutually 
vertex- disjoint paths connecting X with Y in the edge weighted network L is at 
least 1-f. 

Proof. To use the isolation lemma, let the underlying set to consist of all edges in 
the network L. Next, note that a set of mutually vertex-disjoint paths connecting 
X with Y achieving a minimum cost consists of simple paths and thus it can 
be identified with the set of edges on the paths. Let P be the family of all sets 
of mutually vertex-disjoint simple paths connecting X with Y in the network 
L. By the setting of new costs c'(e), solely those sets in P that achieved the 
minimum cost, say D, under the original costs c(e) can achieve a minimum cost 
under the new costs c'(e). So, we can set F to the aforementioned sub- family 
of P, and define the weight of a set of k paths in F as the sum of the weights 
w(e) of the edges e on the paths in this set in order to use the isolation lemma. 
By the isolation lemma, there is a unique set S in F that achieves the minimum 
weight w(S) with the probability at least 1 — — . The corollary follows since each 
set S in F has the cost c'(S) equal to D + w(S). □ 

Throughout the rest of this section, we shall assume that each of the m edges 
e in the network L is assigned the cost c'(e) as in Corollary 1 and that r € [n ^]. 



Suppose that we know the minimum cost of a set of k mutually vertex-disjoint 
paths connecting X with Y in the network L with the edge costs indicated, and 
such a minimum-cost set is unique. Then, it is sufficient to show that we can test 
quickly in parallel if the network L with an arbitrary edge removed still contains 
a set of k mutually vertex-disjoint paths connecting X with Y that achieves the 
minimum cost. By performing the test for each edge of L in parallel, we can 
determine the set of edges forming the unique minimum-cost set of k mutually 
vertex-disjoint paths connecting X with Y. 

To carry out these tests, we need to generalize the polynomial CQl,u to a 
polynomial CPL, e ,u, where e is an edge in L and U is a cost constraint from 
[mr(mC + 1)] = [Cn°^]. Let Hl&u be the family of all proper sets of k walks 
in the network L with the edge e removed that have total at most U. (In the 
total cost of a set of walks, we count the cost of an edge the number of times 
equal to the sum of the multiplicities of the edge in the walks.) 

As in the definition of Q L l assign a distinct variable x e to each edge e in L, 
and for a walk W € H L e ^ w , let M w be the monomial, where x e has multiplicity 
equal to the number of occurrences of e in W. The polynomial CPl^,u is defined 

h yT,seH Lte , w Uwes M w- 

By using the proof method of Lemma 1, we obtain the following counterpart 
of this lemma for CPL, e ,u- 

Lemma 8. For the edge-costed network L with m edges, edge e, and U G 
[Cn°W], there is a proper set of k mutually vertex-disjoint walks of total cost 
< U in the network L with the edge e removed iff CPl.^u is not identical to 
zero over a field of characteristic two. 

Next, we obtain the counterpart of Lemma 3 for CP^^jj following the proof 
of Lemma 6. 

Lemma 9. CPl. £ ,u can be evaluated for a given assignment of values over a 
field F 2 o (log n ) of characteristic two in 0(k log(Cn) +log 2 (Cn)) time by a PRAM 
using 2 k (kCn)°^ processors. 

Now, we are ready to derive our main result in this section. 

Theorem 3. There is a randomized PRAM returning almost certainly (i.e., 
with probability at least 1 — where a > 1) a minimum-cost set of k mutually 
vertex- disjoint paths connecting X with Y in the network L with the original 
edge costs in [C] (iff such a set exists) in 0(k\og(Cn) + log 2 (Cn)) time using 
2 k (kCn)°( 1 ) processors. 

Proof. We set r to, say, n 2 m, and specify the new edge costs c'(e) in the network 
L drawing the weights w(e) uniformly at random from [r] as in Corollary 1. Next, 
for each U <G [mr(mC +1)] = [Cn°W], we proceed in parallel as follows. For each 
edge e of the network L, we test the polynomial CPL. e ,u for the non-identity 
with zero by using Lemma 4 and Lemma 9 (we can perform a linear in n number 
of such tests in parallel in order to decrease the probability of the one-sided error 



to an exponentially small one). Next, we verify if the edges that passed the test 
positively yield a set of k mutually vertex-disjoint paths connecting X with Y. 
For example, it can be done by checking for each cndpoint of the edges outside 
X U Y if it is shared by exactly two of the edges, and then computing and 
examining the transitive closure of the graph induced by the edges (see [19]). If 
so, we save the resulting set of paths of total (new) cost < U. By Corollary 1, 
there is a U € [Cn°^] for which the above procedure will find such a set of 
paths that achieves the minimum (original) cost with probability at least 1 — K 

a 

6 Minimum-cost logarithmic integral flow is in RNC 2 

The following lemma is a straightforward generalization of a folklore reduction 
of maximum integral flow to a corresponding disjoint connecting path problem 
(for instance cf. [7]) to include minimum-cost integral flow. We shall call a flow 
proper, if it ships each flow unit along a simple path from the source to the sink. 

Lemma 10. The problem of whether or not there is a proper integral flow of 
value k and cost D from a distinguished source vertex s to a distinguished sink 
vertex t in a directed network with n vertices, integral edge capacities and edge 
costs in [C] can be (many-one) reduced to that of whether or not there is set of k 
mutually vertex-disjoint simple directed paths of total cost D* , where |_fr~J = D, 
connecting two distinguished sets of k vertices in a directed network on 0(kn 2 ) 
vertices in (9(logfc + logn) time by a CREW PRAM using 0(kn 2 + k 2 n) pro- 
cessors. 

Proof. Let K — (V,E) be the directed network with integral edge capacities, 
edges costs in [C] and the distinguished source vertex s and sink vertex t. Since 
we are interested in a flow of value k, we can assume w.l.o.g that all edge capac- 
ities do not exceed k. 

We form a directed network K* on the basis of the network K as follows. 

Let v G V. Next, let E in (v) be the set of edges in K incoming into v, and 
let E out (v) be the set of edges in K leaving v. For each e e E in (v) and i e 
[capacity (e)], we create the vertex Vi n {e,i). Analogously, for each e' € E out (v) 
and i' € [capacity (e 1 )}, we create the vertex v ou t(e',i'). Furthermore, we direct 
an edge from each vertex Vi n (e, i) to each vertex v out (e',i'). To each such an 
edge, we assign the cost 1. Also, for each edge e = (v, w) of K, we direct an edge 
from v ou t(e,i) to Wi n (e,i) for i e [capacity (e)]. To each such an edge, we assign 
the cost c(e)kn. See Fig. 2. Let X' be the set of vertices of the form s out (...), 
and let Y' denote the set of vertices of the form t j n (...). Create an additional 
set X of k vertices and from each vertex in X direct an edge to each vertex in 
X'. Symmetrically, create another additional set Y of k vertices and from each 
vertex in Y' direct an edge to each vertex in Y. 

It is easy to observe that there is a proper integral flow of value k and cost 
D from s to t in the network K iff there is a set of k mutually vertex-disjoint 




Fig. 2. An example of a vertex of the network K and the corresponding part of the 
network K* . 

simple paths of total cost D* connecting X with Y in the network K* , such that 

L^J =d. 

Now it is sufficient to note that the construction of K* , X and Y on the 
basis of K easily implemented by a CREW PRAM in 0(\ogk + logn)-time 
using 0(kn 2 + k 2 n) processors, where n is the number of vertices in K. □ 

By combining Theorem 1 with Lemma 10, we obtain our first main result. 

Theorem 4. The minimum cost of a flow of value k in a network with n ver- 
tices, a sink and a source, integral edge capacities and positive integral edge 
costs in [C] can be found by a randomized PRAM, with errors of exponentially 
small probability in n, running in 0(k\og(Ckn) + log 2 (Ckn) time and using 
2 k (kCn)°( 1 ) processors. 

By combining in turn Theorem 3 with the finding variant of Lemma 10 using 
exactly the same reduction, we obtain our second main result. 

Theorem 5. There is a randomized PRAM algorithm returning almost cer- 
tainly a minimum-cost flow of value k (iff a flow of value k exists) in a network 
with n vertices, a sink and a source, integral edge capacities and edge costs in 
[n°^], in 0(k\og(kn) + \og 2 {kn)) time using 2 k (kn)°^ processors. 

Corollary 2. The problem of finding a minimum-cost flow of value O(logn) 
in a network with n vertices, a sink and a source, and integral edge capacities 
bounded polynomially in n admits an RNC 2 algorithm. 

7 Final remarks 

We have resented a new approach to the minimum-cost integral flow problem. 
In particular, it yields an RNC 2 algorithm when the flow supply is (at most) 
logarithmic in the size of the network. 



All our results can be extended to include undirected networks by a straight- 
forward reduction. 



References 

1. R.K. Ahuja, T.L. Magnanti, and J.B. Orlin. Network Flows: Theory, Algorithms, 
and Applications. Prentice-Hall, Inc., 1993. 

2. R.J. Anderson and J.C. Setubal. On the Parallel Implementation of Goldberg's 
Maximum Flow Algorithm. Proc. Annual ACM Symposium on Parallel Algorithms 
and Architectures (SPAA), pp. 168-177, 1992. 

3. A. Bjorklund. Determinant sums for undirected Hamiltonicity. Proc. 51th IEEE 
Symposium on Foundations of Computer Science (FOCS), pp. 173-182, 2010. 

4. A. Bjorklund, T. Husfeldt and N. Taslaman. Shortest Cycle Through Specified 
Elements. Proc. Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), 
2012. 

5. R.A. DeMillo and R.J. Lipton. A probabilistic remark on algebraic program testing. 
Information Processing Letters 7, pp. 193-195, 1978. 

6. J. Edmonds. Systems of distinct representatives and linear algebra. J. Res. Nat. 
Bur. Standards Sect. Vol. 71B, No. 4, pp. 241-245, 1967. 

7. S. Even. Graph Algorithms. Computer Science Press Inc. 1979. 

8. J. Edmonds and R.M. Karp. Theoretical improvements in algorithmic efficiency 
for network flow problems. Journal of the ACM 19(2), pp. 248264, 1972. 

9. L.R. Ford jr. and D.R. Fulkerson. Maximal flow through a network. Can. J. Math 
8, pp. 399-404. 

10. D.R. Fulkerson. An out-of-Kilter method for minimal cost flow problems. SIAM 
J. Appl. Math. 9, pp. 18-27. 

11. Z. Galil and V. Pan. Improved processor bounds for algebraic and combinatorial 
problems in RNC. Proc. 26th IEEE Symposium on Foundations of Computer 
Science (FOCS), pp. 490-495, 1985. 

12. A. Goldberg. Parallel Algorithms for Network Flow Problems. In Synthesis of 
Parallel Algorithms, editor J.H. Reif., Morgan-Kauffman, 1993. 

13. L.M. Goldschlager, R.A. Shaw and J. Staples. The Maximum Flow Problem is Log 
Space Complete for P. Theoretical Computer Science 21(1), pp. 105-111, 1982. 

14. R. Karp, E. Upfal, and A. Wigderson. Constructing a perfect matching is in 
random NC. Combinatorica 6(1), pp. 35-48, 1986. 

15. I. Koutis. Faster algebraic algorithms for path and packing problems. Proc. 
35th Annual International Colloquium on Automata, Languages and Program- 
ming (ICALP), Lecture Notes in Computer Science 5555, pp. 653-664, 2009. 

16. E.L. Lawler. Combinatorial optimization: Networks and matroids. Holt, Rinehart 
and Winston, New York, NY. 

17. K. Mulmuley, U.V. Vazirani and V.V. Vazirani. Matching is as Easy as Matrix 
Inversion. Combinatorica 7(1), pp. 105-113, 1987. 

18. J.B. Orlin and C. Stein. Parallel Algorithms for the Assignment and Minimum- 
Cost Flow Problems. Operations Research Letters, 14, pp. 181-186, 1993. 

19. J.H. Reif (editor). Synthesis of Parallel Algorithms. Morgan-Kauffman, 1993. 

20. V.V. Vazirani. Parallel Graph Matching. In Synthesis of Parallel Algorithms, 
editor J.H. Reif., Morgan-Kauffman, 1993. 

21. J.T. Schwartz. Fast probabilistic algorithms for verification of polynomial identi- 
ties. Journal of the ACM 27(4), pp. 701-717, 1980. 

22. R. Williams. Finding paths of length k in 0*(2 k ). Information Processing Letters 
109, pp. 301-338, 2009. 



